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Abstract. In Au+Au collisions at RHIC, elliptic flow was found to rapidly 
decrease as a function of rapidity. We argue that the origin of this phenomenon is 
incomplete thermalization of the initial fireball outside the midrapidity region and 
show that it can be quantitatively related to the analogous effect at midrapidity 
in peripheral collisions and in collisions at lower beam energies. 



1. Introduction and overview 

Semicentral Au+Au collisions at RHIC exhibit very strong elliptic flow w 2 which 
exhausts the prediction from ideal fluid dynamics (for a review see PQ). This has been 
interpreted as strong evidence for early thermalization of the collision fireball, on a 
timescale of < 1 fm/c and at energy densities of ~ 10 — 20 times the critical value for 
quark-gluon plasma formation 0. However, in more peripheral Au+Au collisions at 
RHIC the measured elliptic flow remains increasingly below the hydrodynamic value, 
and in lower energy Pb+Pb and Au+Au collisions at the SPS and AGS it does not 
reach the hydrodynamic limit even in central collisions H] . The STAR and NA49 
Collaborations have suggested (see Fig. 22 in [3] and Fig. 25 in 0]) that the discrepancy 
between the hydrodynamically predicted and measured elliptic flow at midrapidity 
scales with (1/S) dN^/dy where dN^/dy is the measured charged multiplicity density 
at midrapidity and S is the initial transverse overlap area between the colliding nuclei. 
When multiplied with the average transverse mass (m±) per hadron and divided by 
To, this scaling variable is just the Bjorken energy density at proper time To [S]. This 
suggests that the degree of thermalization reached early in the collision (as measured 
by the fraction i;™ eas /i!2 ydro of the hydrodynamically predicted elliptic flow achieved 
in the experiment) is controlled by the initial energy or particle density. This makes 
sense since the collision rate is proportional to the density of scatterers. 

The PHOBOS Collaboration first observed jHj that V2 shows a strong pseudo- 
rapidity dependence, with a shape that seems to roughly follow the charged particle 
pseudorapidity distribution. This has now been confirmed by STAR 0. BRAHMS 
and PHOBOS also reported that the pseudorapidity distributions of both iV c h [HI OH 
and t>2 EH have very similar shapes for all collision centralities. The similarity in 
shape between the rapidity distributions of V2 and iV c h becomes even stronger if the 
Jacobian between rapidity y and pseudorapidity rj (which affects V2(r)) and dN^/dr] 
mostly near midrapidity and with opposite signs is taken out. 

Hydrodynamic models can reproduce the shape of dN^/dr] but not that of V2(t]) 
|12j . Instead of falling off with increasing rapidity, the hydrodynamically predicted 
1)2(11) actually first increases with |?7| before eventually dropping steeply near the 
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projectile and beam rapidities ^3 EH- ^ e ^ n( ^ * w0 reasons f° r this behaviour: (i) 
The initial transverse spatial deformation e x of the nuclear overlap region, which 
causes the anisotropic pressure gradients driving the buildup of elliptic flow, slightly 
increases with space-time rapidity \r] s \ = \ In favoring larger v-i at 77 7^ 0. (ii) The 
initial transversally averaged energy density e(ri s ;T ) = (e(r, ?] s ; r )) decreases with 
\r) s \. This reduces the softening effects on the flow buildup from the quark- hadron 
phase transition and again results in larger V2- The same effect leads to larger V2 at 
midrapidity when one reduces the beam energy from RHIC to SPS and AGS energies 
|14| . V2 falls to zero only near the projectile and target rapidities since there the initial 
energy density becomes so small that no time is left for flow buildup before freeze-out. 

Hirano interpreted the observation that for r\ ^ the measured 1*2 in Au+Au at 
RHIC is significantly below the hydrodynamically predicted one as evidence that local 
thermal equilibrium in the initial state of the collision is only achieved at midrapidity 
|T2). We here support this conclusion by showing that the ratio i;™ eas /i>2 ydro scales 
with the initial transversally averaged energy density as a function of space-time resp. 
pseudo-rapidity in the same way as previously observed by STAR and NA49 as a 
function of centrality at midrapidity E]- Both effects indicate that the initial 
thermalization of the system and the validity of hydrodynamics are largely controlled 
by the initial energy or particle density achieved in the collision. 



2. Initialization and longitudinal structure of the reaction zone 

Understanding the rapidity dependence of elliptic flow requires a model which 
correlates the initial transverse distribution with longitudinal position, i.e. with space- 
time rapidity r/ s . We use Bjorken scaling r) s = y [S] during the particle formation stage 
to relate initial longitudinal position to initial longitudinal momentum by z = v^t on 
a surface tq = (i 2 — z 2 ) 1 / 2 = const. The average longitudinal momentum is correlated 
with transverse position r by overlap geometry coupled with momentum conservation 
|15| : At transverse position r two cylinders of matter with thicknesses ti t 2(f;b) 
= TAu(r±b/2)% hit each other, carrying total energy 7beam(ii+i2) and net z-momen- 
tum 7beami'beam(ii^^2)- In a collision with impact parameter b, the produced matter 
thus has a rapidity distribution centered at the r-dependent avarage rapidity |15j 

w ^ 1 ^ (ti(r;b)+t 2 (r;b)) + 

^bcam 

(t 1 (r;b)-t 2 (r;b)) ... 
{Vl ' 2 1 \t 1 (v-b) + t2(r-b))~v hc&m {t 1 {r-,b)-t 2 {r-b))- {) 

In coordinate space the center of mass of the matter produced at transverse position 
r is located at rjf^fr) = Y(r). We distribute the initial matter symmetrically around 
this mean space-time rapidity as follows: Noting that 1+1 and 3+1 dimensional 
hydrodynamical calculations |161 117) show very little dynamical evolution in rapidity 
direction, we simply identify the shape of the initial entropy distribution dS/dn s 
with that of the measured final charged particle rapidity distribution dN^/dy. We 
parametrize the latter by a double Gaussian, 

(v-») 2 («+°) 2 . , 

dN ch /dy a. e + e , (2) 

with a free width parameter a. Combining Eqs. 11121) and noting that the normalization 
of the rapidity densities at midrapidity scales roughly with the number of wounded 

7a (r) = J°° dz pa (r, z) is the thickness function of a nucleus of mass A, calculated from a Woods- 
Saxon profile PaM = Po/[ ex p(( r— Ra)/&) + l]j with i?Au = 6.37fm and <5 = 0.54fm. 
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Figure 1. Left: Pseudorapidity spectra from 200 A GeV Au+Au collisions with 
different impact parameters. PHOBOS data |5| are compared with a double 
Gaussian parametrization (see text). Right: Hydrodynamic time evolution of the 
momentum anisotropy e p at space-time rapidities r] s = 0, 1, 2, 3, 4. 

nucleons, we arrive at the following form for the initial entropy distribution: 

s(r, b) oc („r (r; b) + (r; 6)) ^- ^^M! + e - ^£^ ) . (3) 

Integration over r gives the rapidity distribution dS/drj s oc dN c h/dy, which can 
be transformed to pseudorapidity, dN c h/drj, by multiplying with the Jacobian 
(p_i_) cosh?7/(m^ + (p^) 2 cosh 2 /;) 1 / 2 , where we treat all charged hadrons as pions with 
(p_i_) w 0.4 GeV. A single overall normalization constant and the width a of the double 
Gaussian are fitted to the data. The left panel in Fig. ^ shows such a fit for 
200 A GeV Au+Au data with a = 1.9±0.05. Similar fits work well for Au+Au data at 
Vs = 130 and 19.6AGeV, with 0= 1.8+0.05 and 1.15+0.05, respectively. Note that 
the double Gaussian ansatz does not work for d+Au collisions. 

3. Time evolution and anisotropic flow observables 

For each value of rj s , we evolve these initial conditions with a locally boost-invariant 
hydrodynamic code ^1] with the equation state of 18 to simulate the transverse 
evolution of the reaction zone. Here we represent minimum bias Au+Au data by a 
fixed average impact parameter of b = 6.8 fm. The evolution is started at To = 0.6 fm/c 
and stopped at decoupling energy density edcc = 0.075 GeV/fm 3 , both independent of 
rapidity. As noted by Hirano |17j , the smaller transverse size of the fireball at forward 
rapidities should lead to an earlier decoupling at higher temperature and smaller 
transverse flow, causing steeper p_L spectra and smaller pj_-integrated elliptic flow at 
forward rapidities. The same point was made previously by Teaney for the smaller 
collision systems formed in peripheral collisions at midrapidity [T^]. This effect is, 
however, not enough to explain the strong rapidity dependence of the measured vi 
|17| . and it is not included in our present calculations. 

In the right panel of Fig. ^ we show the time evolution of the momentum anisotro- 

IrpXX rpyy \ 

py £p = ^T^^+Tyy) a * var i° us rapidities ?7 S .|| As n s increases, the initial energy density 

I Beyond r — to — 11 — 12 fm/c our numerical results for e p are affected by matter leaving our finite 
grid area and tend to be too large. From 1141 we know that e p saturates beyond r — to ~ 12 fm/c. 
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Figure 2. Left; Differential directed flow Ui(px) of directly emitted pions (no 
resonance decays) for rj s = y = 0, 1, 2, 3, 4. Except for a region of positive v\ at 
< p± < 0.5 GcV and a shift of the rest of the curves by about 0.5 GeV to larger 
pj_, the curves for direct protons look similar. Right: p^-integrated elliptic flow 
V2 for direct pions as a function of pseudorapidity r], compared with data for all 
charged particles [5][7|. The hydrodynamic 1)2 values have been corrected with an 
energy density dependent "thermalizations coefficient" as described in the text. 
The Jacobian for the transformation from y to rj has been included. 



at r] s decreases, and the time evolution of e p follows the same pattern as previously 
observed at midrapidity when reducing the collision energy (see Fig. 7 in 14 ). 

At forward rapidities the transverse overlap region becomes asymmetric and is 
shifted sidewards in the x (or impact parameter) direction. This turns out to give 
rise to a non-zero directed flow signal vi(p±) which increases with |ry s | (left panel in 
Fig. |2J. Of course, since the colliding matter receives no overall transverse kick, the 
px-integrated directed flow is zero. 

The hydrodynamically calculated elliptic flow v 2 {rj) has the same general shape as 
previously obtained by Hirano with a fully (3+l)-dimensional code. We now correct 
this hydrodynamic behaviour with a "thermalization coefficient" F(x) which is fitted 
to midrapidity data in peripheral and lower-energy collisions [3] 0] . F depends on 
the initial transversally averaged energy density at rapidity y = rj s through the ratio 
x(rj s ) = (e(ry s ))/eo (where eo = 9.5 GeV/fm 3 is the average initial energy density in 
central Au+Au collisions at 130^4 GeV). As discussed in the Introduction, this scaling 
variable is, up to a multiplicative constant, identical with the variable {\/S)dN/dy 
found by STAR and NA49 to control the magnitude of v 2 at midrapidity 0] . We 
parametrize the behavior shown in Fig. 25 of with a simple linear function F(x) = 
^ y dro = 0.15 + 0.85 x for x < 1 while F(x) = 1 for x > 1. (x = 1 corresponds in Fig. 25 

V 2 

of [U to (l/S)dN ch /dy = 25fm- 2 .) The corrected vf eas (rj) =F(x(t])) ■ v% ydm (r)) 
for 6 = 6.8fm is shown by the full circles in the right panel of Fig. together with 
minimum bias data from PHOBOS and STAR. Even if our v 2 values are still a bit high 
at | rj | > 2, we see good qualitative agreement with the data. We conclude that the 
same incomplete thermalization effects previously seen at midrapidity in peripheral 
and lower-energy collisions also describe qualitatively the rapid decrease of v 2 at non- 
zero rapidity in minimum bias collisions at RHIC. Local thermalization seems to be 
driven by the local initial energy density reached in the collision. 



Rapidity dependent momentum anisotropy at RHIC 



5 



4. Conclusions and further tests 



Our analysis suggests that the critical initial energy density, which is required for 
reaching sufficiently complete local equilibrium during the early collision stages to 
validate a hydrodynamic description, has for the first time been reached at midrapidity 
in almost central Au+Au collisions at top RHIC energy. At even higher collision 
energies the elliptic flow is therefore expected to follow the hydrodynamic predictions 
over a wider range of rapidities and centralities. However, there is no need to wait for 
the LHC to verify this prediction. One way to test this already at RHIC energies would 
be to explore fully central collisions ( "zero spectators" ) between Uranium nuclei which 
provide slightly larger energy densities than central Au+Au collisions and a significant 
transverse deformation and elliptic flow signal even at zero impact parameter, due 
to the deformation of the projectiles flj. Such full-overlap U+U collisions yield 
very large, transversally deformed fireballs which would also provide a testbed for 
studying more efficiently the predicted nonlinear pathlength dependence of jet energy 
loss (by studying jets as a function of their emission angle relative to the reaction 
plane as extracted from v^) and exploring this energy loss out to larger px values than 
presently possible with the much smaller deformed fireballs created in non-central 
Au+Au collisions. 
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